Laser-induced collective excitations in a two-component Fermi gas 
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We consider the linear density response of a two-component (superfluid) Fermi gas of atoms when 
the perturbation is caused by laser light. We show that various types of laser excitation schemes can 
be transformed into linear density perturbations, however, a Bragg spectroscopy scheme is needed 
for transferring energy and momentum into a collective mode. This makes other types of laser 
probing schemes insensitive for collective excitations and therefore well suited for the detection of 
the superfluid order parameter. We show that for the special case when laser light is coupled between 
the two components of the Fermi gas, density response is always absent in a homogeneous system. 
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I. INTRODUCTION 



Fermionic atoms have already been cooled well below 
the Fermi temperature . Once the quantum regime has 
been achieved, the predicted BCS superfluid transition 
H for two different hyperfine states with attractive inter- 
actions between them would be the next goal to achieve. 
There are several methods proposed for probing the su- 
perconducting gap |^-^] . Some of them are based on the 
use of laser light. For instance, it has been proposed to 
measure the superfluid coherence and the existence of the 
superfluid gap by the absorption of almost on-resonance 
laser light coupled between the two hyperfine states Q or 
between one of them and another atomic hyperfine state 
H . Also the use of off-resonant (or very low intensity) 
laser light Q has been proposed. There both hyperfine 
components are coupled to other excited states not in- 
volved in the pairing, and the excited state populations 
remain negligible. 

In this paper we review and analyze the linear density 
response of the system to the different laser excitations in 
the collisionless regime where the lifetime of the quasipar- 
ticles is much longer than the period of the applied field. 
The motivation for studying the response is two-fold. 
First, in the above-mentioned probing schemes for the 
gap a below-gap density response (e.g. Anderson- Bogoli- 
ubov phonon) would be undesireable because the schemes 
rely on the absence of below-gap excitations. Therefore 
it is of interest to clarify in what kind of laser probing 
schemes density reponses can be avoided. Second, the 
density response itself may be the object of study and 
the laser excitation schemes that can produce a response 
are of interest. 

Off-resonant laser excitation caused by two intersect- 
ing waves which form a density grating in space and time 
(Bragg spectroscopy) has been used to induce collective 
excitations in atomic Bose-Einstein condenstates 0. In 
case of near-resonant laser excitations, we transform the 
non-diagonal atom-light interaction Hamiltonian into a 
form of a density perturbation. We show that the initial 
Hamiltonian of the interacting two-component Fermi gas 



is preserved by this transformation in case of a contact 
interaction and therefore a standard linear density re- 
sponse calculation for a BCS-system can be applied. In 
general, an intensity grating in space and time is needed 
for providing momentum and energy transfer for den- 
sity perturbations. We show that when the laser light is 
coupled between the two paired hyperfine components, 
phonons cannot be excited even when the light forms an 
intensity grating. 

In section II, the Hamiltonian for a two-component 
Fermi gas is given and the density response is defined. 
The laser probing schemes considered in this article are 
summarized in section III. In section IV we transform 
the perturbation caused by the probing lasers into a den- 
sity perturbation by appropiate rotation of the coupled 
states. The results are presented in detail in section V 
and summarized in section VI. The linear response cal- 
culation is presented in the Appendix. 



II. THE SYSTEM 

The Hamiltonian for the two-component Fermi gas, in- 
cluding only 2-body interactions between the atoms reads 



H = J dr J2 



V>£(r) 



2m 



+ V a (r) - \i a ip a (r) 



+ \ E l dr /rfr / 3a/3(r-r')^(r)^(r')V'/3(rOV'a(r) ) (1) 

where the two hyperfine states trapped are denoted by f 
and {. When the interaction gn(r — r') is attractive, the 
system has been predicted to undergo a BCS transition 
|U associated with the appearance of an order parameter. 
We assume the system to be below the critical temper- 
ature and the Hamiltonian H treated within the BCS 
approximation. 

The density response of the system to perturbations, 
caused by e.g. probe lasers or the trapping fields, is stud- 
ied in the linear response regime j^] . The density pertur- 
bation is a linear function of an external potential U(r, t) 
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expressed through the density-density response function 
xM;r',i'): 



Sn(r, t) = 5n-\ (r, t) + 8n[ (r, t) 



(2) 



where the density operators are defined in the usual way 
n Q= | j |(r) = ■03 1; (r)?/; Q (r). Fourier transforming the previ- 
ous expression one obtains 



«»(k,w) = x(k,w)^(k,w) 



(3) 



The poles of the reponse function x give the collective 
modes of the system while its modulus gives the spectral 
weight of the modes. 

It is known that a perturbation of the form 
U(t, t) [V>j (r)ipi (r) + ip^(r)ipi(r)] leads in the long wave- 
length limit to appearance of the Anderson-Bogoliubov 
phonon § (fl], that is, a collective mode with energy be- 
low the superconducting gap. 

In this paper we study the reponse function x(k, 
and potential t/(k, lo) for a perturbation H' created by 
different laser probing schemes. Most of the schemes have 
been considered in the literature as proposed techniques 
for observing the superconducting gap. 



III. PROBING SCHEMES 

A. Laser coupling between two states 

a) Coupling the two paired states. Let us assume that 
the states f and j are coupled by ligh t fl . Using the 
Rotating- Wave- Approximation (RWA) 112] ", the interac- 
tion of the laser light with the matter fields can be de- 
scribed by a time-independent Hamiltonian in which the 
detuning 8 plays the role of a chemical potential and 
a Rabi frequency fl(r) characterises the local strength 
of the interaction. The perturbation Hamiltonian is 
H' = + Ht, where the transfer Hamiltonian (Ht) 
and are given by 



drf)(r)V}(r)Va(r) + fi*(r)^J(r)^ T (r) 



II, = I [ drV4(r)V T (r) - | / dv^Av^^v) (4) 



In principle, one can imagine a direct coupling of the 
states | and J, by one field, in which case S and fl cor- 
respond directly to that field. In practice, for coupling 
of the two hyperfine states which have an energy dif- 
ference far from any laser frequency one may prefer a 
Raman transition. In that case the Rabi frequencies and 
6 depend on the parameters from both lasers, especially 

fi cx n$n 2 - 

b) Coupling to a non-paired state. When one of the 
states a — {f or J,} is coupled to some excited hyperfine 
state e H , the perturbation reads 



//,, = I [ dr^(r)Va(r) - £ / dr^( r )^ e ( r ). ( 5 ) 



c) Far off-resonant light. In H|J the coupled light was 
proposed to be almost on-resonant, that is, atomic pop- 
ulation is transfered in the excitation process. In || cou- 
pling of the paired states with some excited states was as- 
sumed to be done by far off-resonant or very weak inten- 
sity light meaning the excited state population remains 
negligible and the excited state can be eliminated from 
the problem. Far off-resonance coupling can be treated 
by starting with an initial Hamiltonian of the same form 
as (||), then after adiabatically eliminating the excited 



states one arrives at —g-^a^a, c.f. section IV 



B. Space and time dependence of the perturbation 



d) Spatial variation of the intensity. The common way 
of creating intensity gratings used e.g. for creating col- 
lective excitations in Bose-Einstein condensates [[7j is to 
use two intersecting waves. When two intersecting waves 
with the same polarization but different wave- vector and 
frequency are coupled to the same two-level system, the 
total Rabi frequency is fl = fl ie ^ r - wlt + fl 2 e lii2r - luJ2t . 
After making the RWA for the frequency Ui, fl = 
fiie* 1 ' 1 " + f7 2 e ik2 - r+i " 12 * where u X2 =u) X - u 2 . 

Also the beam profile may vary spatially in the cloud 
size scale. The Rabi frequency would then be for instance 
fl(r) cx e~ r /( 2a ' for a Gaussian beam shape. 

As will be shown in the following section, the essen- 
tial feature considering a possible density response is 
whether |f2| 2 is time and space dependent. Clearly, for 
Bragg spectroscopy as considered above this is true. For 
the Gaussian beam profile, |f2| 2 is dependent on posi- 
tion but not on time. In contrast, for a single laser 
fl = |S!|e 4k ' r-lw * and for a Raman excitation fl oc 
\fl 1 \\fl 2 \e ik ^ r - lullt e-^- r+liA)2t , thus in both cases \fl\ 2 is 
a constant. 

e) Time dependece of the perturbation. As discussed 
above, Bragg scattering schemes lead to a time- varying 
intensity grating, that is, the energy coi 2 is transferred to 
the system via a perturbation that is proportional to the 
intensity. 

Another possible source of time-dependence is the 
turning on of the perturbation. As is usually done in 
the linear response theory we consider smooth turn- 
ing on of the perturbing potential. 

Direct intensity modulation of the probing laser(s) 
would be one more source of time-dependence. 
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IV. TRANSFORMING H T INTO A DENSITY 
PERTURBATION 

In order to use linear response theory we have to con- 
vert all the perturbing Hamiltonians of the previous sec- 
tion into a density perturbation. 

a) The perturbation term Ht involves products of the 
field operators of the different atomic states. In the 
mean field approach they correspond to the Fock terms 
< i^I(t)iPi(t) > which are zero. Therefore no linear re- 
sponse theory can be directly applied. 

A perturbation of the type (note that the two- 
component ip here is not the same as that in the standard 
BCS-theory, used for instance in the Appendix) 



^{v)Wi){v) 



5/2 
fi*(r) 



n(r) 

-5/2 



(0) 



can be digonalized for any 5 and f2 by an appropiate ro- 
tation in the space of the two states j and J,, ip = Uip 
and U=R z {a)R v {[3)R z {"i) with Rj(d) = exp(io-j6), j = 
{x, y, z} and o~j is the corresponding Pauli spin matrix. 
The transformation matrix is 



U 



e i(«+7)/2 C0S/ 3/2 e t{a -^/ 2 s : m(3/2 \ 

_ e - 4 (a- 7 )/2 sin/3 / 2 e -«("+7)/2 cos/ 3/ 2 ) 



(7) 



The rotated perturbation matrix is U^WU. The off- 
diagonal terms are 

e-^((S/2) sin(/3) - e la il* sin 2 (/3/2) + Q cos 2 (/3/2) e - m ) 

(8) 

and the complex conjugate. Off-diagonal terms that are 
zero can be obtained by a rotation with the Euler angles 
{7 = 0,cosa = %^,tan/3 = -^} which yields 



(I) 2 



nn* 











+ (l) 2 



(9) 



for the perturbation matrix. This procedure is exten- 
sively used in describing the interaction of laser light with 
a two-level system. The form |0| 2 /<5 for far off-resonant 
light comes from (§) for \5\ >> 

One has to check whether the rotated states feel 
the same two-body interaction that the non-rotated 
ones, that is, whether the rotated states are still de- 
scribed by a Hamiltonian of the initial form which im- 
plies BCS pairing. Assuming, for simplicity, that ijj-\ 
and ipi see the same potential V a — fi a , we obtain 
V4(r)<M r ) + V>i(r)^(r) i>\(r)i> r (r) + ^j(r)^(r). 
The kinetic energy term transforms as V^|(r)V'i/'-f (r) + 

V^WVV^W -> V^}(r)V^ T (r)+V^j(r)V^(r) if the 
momentum k transfered by the laser is small. The extra 



terms coming in the transformation of the kinetic energy 
are ~ k,k 2 . The momentum of the laser and the recoil 
energy are very small compared to the momentum of the 
atoms participating in the process (typically close to the 
Fermi-surface, ~ kp). 

The interaction term in the Hamiltonian transforms 
■0| (r) (r')V'j. (r')-0^ (r) — > ^j^j^^f in case of a contact 

interaction g(r — r') cx 5(r— r'). As ip a is a linear combi- 
nation of the two species t and |, one might expect inter- 
action terms between the same species for the new mat- 
ter fields (interactions between both ip] and tp^, and be- 
tween tp^ and tpi). But the interactions of the type ^) a ^> a 
are forbidden by the fermionic behaviour (commutation 
rules are preserved by the rotation) . The requirement of 
a contact interaction is obvious here: tpa(r)tlj a (r) = 0, 
but ip a (r)ipa(r') can be non-zero. As a summary, since 
the fermionic field commutator relations are preserved by 
the rotation, for contact interaction Hq = Hq. 

We have thus shown that the system can be trans- 
formed in such a way that one can study the collective 
mode spectra by linear response theory for a BCS sys- 
tem for a perturbation of the type U(n-\ — njj, where 



U(k,u) = T 



^MMr,i)* + (|) 5 



and the rest of 



the Hamiltonian still corresponds to the BCS Hamilto- 
nian. The main difference to the standard linear density 
response treatment for an BCS system is the minus sign 
between rif and ny. We will show that this leads to zero 
overall response in the homogeneous case. 

b) When diagonalizing the interaction Hamiltonian of 
Eq.(||) in the same way as in a), the initial Hamilto- 
nian Hq is not preserved. We thus obtain a perturbation 
U(n a — n e ) with t/(k, uj) as given above, but with the 
initial Hamiltonian modified, Hq =/= Hq. 

c) In the case of coupling to an excited state by far- 
detuned light the rotation does not modify the initial 
Hamiltonian because at the limit \5\ >> |f2| the trans- 
formed states are very close to the initial ones, ip a ~ ip a 
and ip e ~ ip e , moreover the excited state population is 
assumed to be negligible. Thus the light-atom interac- 
tion Hamiltonian is in the form of a density perturbation, 
and the initial Hamiltonian is preserved. 

In all these cases, U (k, uj) can be either a constant 
U(k, uj) = U' or dependent on k and uj. The latter hap- 
pens when the intensity |fi| 2 is space and time dependent. 
Note that the k-dependence r2e lkr is not sufficient, that 
is, the laser momentum given in a single-laser or Raman 
two-level excitation (near-resonant or off-resonant) is not 
enough to transfer momentum through a density pertur- 
bation. 



V. RESULTS 
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A. GRPA for the case of light coupled between the 
two paired states 

The response of the system to the perturbation 
U(k,uj)(n^ — njj in a homogenous system is calculated 
in the Appendix following the derivation in Jlp| l . The re- 
sult is that for the homogeneous case, the responses of 
two components are equal but with opposite signs and 
no collective modes are excited at any frequencies. 

The response of one of the spin components is L%i = 

Ti-^tT-^r^ri^ i where we use the integrals a, b, c, d, R 
(i-goR)(go+a+bV^ o 

as defined in [hfl. The reponse of the spin component 
has the factor (a+b) which at T=0 in the low q,w region 
is a + b sa -N(e F )c 2 B q 2 / (12 A 2 ) where N(e F ) is the den- 
sity of states at the Fermi surface and the sound velocity 
c? B = v F /3. This factor is very small and it makes the 
response for one spin component negligible even without 
the cancellation effect. This is because the superfluid 
cannot participate in the relative motion of the two spin 
components 14 1. Above T c , b vanishes and the response 
of one of the spins is not negligible anymore, as the in- 
tegral a is some fraction of N(€f)- However, due to the 
cancellation the overall response is zero. 

We have considered the homogeneous case, that is, 
when the trapping potential treated within the local den- 
sity approximation is sufficient. The strongly trapped 
case can be approached by using multipole expansions 
|^5| : the overall response is not zero for spin-dipole exci- 
tations in the trapped case above T c . The difference to 
our results is a consequence of the geometry, as in the 
homogeneous case we get the same spatially uniform re- 
sponse in opposite directions for both spin species. As 
T — > 0, one gets vanishing response also for the inho- 
mogeneous case Q , at least at energies below the gap, 
because one needs to break pairs in order to have relative 
motion of the two components. 



B. Density response for different probing schemes 

a) Coupling the paired states. The perturbation reads 
U(n^ — rijj for the rotated states, where U(k, uj) — 
^F[^/|0(r, t)\ 2 + 8 2 ]. When the light is far off-resonant, 
one can take the limit U(h,u>) ~ .F[|f2(r, t)\ 2 /S\. In 
the case of homogeneous laser intensity ft = |£l|e lk ' r , 
U cx <J(aj)£(k). There is no density reponse simply be- 
cause U does not provide the momentum and energy for 
a collective excitation. 

A different behaviour arises when using a Bragg scat- 
tering scheme to provide an intensity grating. Assum- 
ing |fix| = (S^! , ri 2 as introduced in section III d)) 
one obtains |f2| 2 = |f2i| 2 (2 + 2cos(wi 2 i — k i2 • r)) ~ 
|fii | 2 2 cos((jJi 2 i— ki2 • r) and the density response x(w, k) 
has to be analyzed at x( w i2>k 12 ). The density response 
X gives zero in the homogeneous case for all tempera- 
tures, as one component annihilates the response of the 



other, and for the inhomogeneous case well below T c |T|] 
as considered in the previous subsection. 

For a general case where the beam profile of the 
laser light gives a spatially varying intensity of the form 

r 2 

|£l(r)| 2 ~ e^, the far off-resonant perturbation reads 

?7(k, u)) = 1 g 1 ' ~ e~^~ . If cr is some fraction / of the 
trap size Rtf, one can estimate the momentum given 
by P — ■ A collective mode with an energy of the or- 
der of the gap energy (e.g. half the gap) has momentum 
q ~ V3^f |,| and the ratio between them is \ = j^. 
For typical parameters, this could be a non- negligible 
number, however, this alone is not enough to produce a 
density response because of the lack of the time depen- 
dence, i.e. £/(w,k) tx <5(w)/(k). 

b) Coupling to an excited state. When diagonalizing 
the interaction Hamiltonian in this case the initial Hamil- 
tonian was modified. Therefore the linear response calcu- 
lation for a BCS system cannot be directly applied. One 
may guess, however, that the arguments about the per- 
turbation U being time- and space-dependent will hold 
also in this case even when the response function \ would 
have a different form. Therefore we do not expect a den- 
sity response except when a Bragg spectroscopy scheme 
is used. 

c) Far off-resonant coupling. When both paired states 
are coupled to some excited states and the light is far de- 
tuned, the density perturbation potential is proportional 
to |f2| 2 /(5. No collective mode of the system is excited if 
there is no intensity variation in space and time. When 
considering intensity modulations of the type |f2| 2 ~ 
|fii | 2 2 cos(cji — k-_r), an Anderson-Bogoliubov phonon 
can be excited [BJlll . 



VI. CONCLUSIONS 

We have reviewed the density response of some typical 
laser excitation schemes. We have shown that since the 
BCS Hamiltonian for a contact interaction is preserved 
under a rotation, most of the considered laser excitations 
can be expressed in terms of a perturbation acting on 
the density. In this form, the perturbation potential is 
proportional to |f2| 2 where f2 is the (effective) Rabi fre- 
quency. Therefore, even when the laser light provides mo- 
mentum and energy (SI cx e lkL ' r e l " I - t ), the transformed 
potential acting as a density perturbation is not time- 
and space-dependent. This leads to absence of a density 
response whenever |f2| 2 is a constant spatially and tempo- 
rally. This makes many proposed laser-probing schemes 
well suited for observing the superconducting gap since 
they do not induce below-gap collective excitations. 

For Bragg scattering, |f2| 2 cx f(r,t). In this case 
Anderson-Bogoliubov phonons can be excited, in gen- 
eral. The exception is the case when the laser(s) cou- 
ple between the two paired components of the gas. We 
have shown that, in the homogeneous case, the den- 
sity response becomes zero because the contributions 
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of the components cancel each other. In a harmonic 
trap, spin-dipole response is predicted for temperatures 
near T c . Therefore, the presence or absence of low- 
energy collective excitations under a perturbation of the 
type U(u>, fc)(ri| — rijj could be used to observe whether 
the trapped system can be approximated by a homoge- 
neous system (local density approximation) or whether 
the trapping effects are dominant. 



VII. APPENDIX 

The formalism used is based on the functional dif- 
ferentiation technique as described in jl3| following the 
derivation in [fTof . The response function \ is viewed as 
a functional derivative of the one-particle matrix Green's 
functions G with respect to the external field U(r, t). To 
allow for pairing in the superconducting phase, it is con- 
venient to work with a single-particle Green's function G 
given by a 2 x 2 matrix defined as 



our case L°(l,2,5) = G(1,5)G(5,2) because the pertur- 
bation matrix W(l) — U(1)X due to the minus sign in 
U(n-\ — rii). ^From now on, we denote by the subindex 
AB the quantities that are calculated for the perturba- 
tion matrix W(l) = U(l)a 3 . 

It is useful to rewrite the GRPA integral equations 
in terms of irreducible two-particle Green's functions 
Zjj(l,2, 5), and for the homogeneous system, to Fourier 
transform in order to solve coupled equations. The result 
reads M 



Lij(q,ioj n ) = Lijiqjojr, 



Labi] (q, iu n )gn (q)£» (q, i^n 

(12) 



where equal indices mean summation over the possi- 
ble values. We denote by the subindex AB the three- 
point correlation functions for the Anderson-Bogoliubov 
phonon type of perturbation. 

The sum of diagonal terms reduces to 



G(l,2) 



-<T*(l)*t( 2 )> = (^l^j 



F(l,2) 
-G x (2,l) 



where "J = (ip^ ^\) T an< i 1 = ( r ii r i)- Imaginary times 
(Matsubara formalism) are used so that one can deal with 
finite temperatures. In the absence of external fields, 
the equal time (j2 — t^) single-particle Green's function 
components reduce to F(r) = (ipi(r)ip^(r)}, the s-wave 
order parameter, and Gj.j(r) = (nj j(r)). 

^From the equation of motion of the Green's func- 
tion one gets the generalized Dyson equation for G(l,2) 
in terms of Go, the non- interacting single-particle 
Green's function, the matrix self-energy S(l, 2) which is 
evaluated in the pairing approximation (Hartree-Fock- 
Bogoliubov), and W(l) which is the the external per- 
turbing field matrix. 

The density response matrix is obtained in RPA by 
taking the functional derivative of the Green's function 
with respect to the external field U. One can define the 



(10) X(q> iuJ n) = £n(q, iu n ) + L 22 (q, iw n ) 

_ Lii{(\, iuj n ) 

1 - 9m{ol)Labii{<\, iw n )' 



(13) 



<5G(1,2) 
SU(5) 



three-point correlation function L(l,2,5) = — 03 

whose limit L(l,2) = L(l,l + ,2) will give the density- 
density response function x(l, 2) = £11 (1, 2) + £22(1, 2). 
Here 03 is the third Pauli matrix. 

When deriving the Dyson equation for G(l,2) with re- 
spect to U in order to get the three-point correlation func- 
tion L, both the self energy matrix and the external field 
matrix W contribute. The lowest-order (single-bubble) 



This shows that all linear perturbations have the same 
poles as the Anderson-Bogoliubov phonon (plus possibly 
some additional ones), but the spectral weight depends 
on the trace of the irreducible two-particle Green's func- 
tions Ljj(q, iijJ n ) of the specific perturbation. 

For a contact interaction (<?-fj.(q) = go), the equation 
for the irreducible correlation function |10[ ] reduces to a 
set of linear algebraic equations 

Lij(q,iuj n ) = Lfjfaiwn) - L ABlklj (q,iuj n )goL M (q,iu; n ), 

(14) 

where we have defined the four-index tensor L ABikl ^ , in- 
dicating components of the two factor matrices. Defining 
column vectors L — ( L\\ L 12 L21 -^22 ) and a 4x4 ma- 



trix L® AB ^ mn as in [|10j], Eq.(14) reduces to 



L=[I + g L° AB }- 1 L° 



(15) 



result 10 for L is given by 



L°(l,2,5) = -03 / d3 / d4G(l,3) 



, 6W{3, 4) 
' 6U(5) 



G(4,2) 



(11) 



For the Anderson-Bogoliubov phonon (see P,10l) the per- 
tubation matrix W(l) — U(l)a 3 , and L^b(1,2, 5) = 
G(l, 5)G(5, 2), where we have introduced G = 03G. In 



As discussed in [[Io| one can get L\- kl by Matsubara fre- 
quency summations, (Eq.(4.27) in JTo[| ) and using the 
symmetry properties then reduce just to 6 independent 
elements. For the weak-coupling limit and contact in- 
teraction, the independent elements are reduced to four: 
a, b, c, d which are integrals defined in jlQ] (Eq.A8). 
lP m can be calculated from the 4x4 matrix by L^T = 
( T° T° T° T° \ 

\ n mi ^iu2 n 2in ^2112 )■ 
Now we derive the 4x4 matrix for our case when 

W^(l) = UT by reconsidering the symmetry properties 

of the matrix elements in order to derive the vector 

that we insert in Eq.(p"5|). The sum of the first and last 

component of L m contribute to Eq. (fl3) and characterize 
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the spectral weight and perhaps some additional poles of 
the response function. Such a perturbation gives from 
(pd|) the lowest order correlation function L°(l, 2, 5) = 
G(1,5)G(5,2) (c.f. £%(1,2,5) = G(l, 5)G(5, 2)). We 
calculate the 4x4 matrix lP mn using the symmetry prop- 
erties of the new L^ kl elements, in the same fashion as 
in jio) (see appendix A there) obtaining 



[1] 



It leads to L°J = 
L° m by [I + g L 





(16) 



[a 

o 

ABmnl 
a+b 



-a - b ) . Multiplying 
as in ( [ilil) one gets i^ = 

This means in = — L22, and 
0. 



^ go+a+b go+a+b 

the spectral weight in (R3h vanishes and \ 

Inserting the well-known result ]lQ,p[ Labu — , 
into Eq.(||), one obtains the spin density response in = 



27V 



go/?. 



(a+i»)(l+ go .R) 



-, where R = A 



(l-g R)(go+a+b)- 

A = a — b are defined in 111 C 



4c z 
1+90 B 



In the case of a general perturbation 



W= U 



1 

c u 



and B = b + d, 



(17) 



where c u is any real number, the single-bubble result is 



1 

c u 



G 

ks / si 



(18) 



By calculating the symmetry properties of 
the elements of L mn , one obtains i^f = 
( a + c u b c — c u c c — c u c —c u a — b ) which leads to 

(c u + l)(a + &) 



in 



^22 — 



and 



ill + i22 — 



l+go(a + b) 
{l-c u )R 



1 + g R 

The reponse for one spin would be 
in 



(1 + g R)(l + g B) 



(19) 



(20) 



(21) 



where 



goAB 



2 (l-c u )+2go(q+h) 
l+g (a+b) 



(a+c u b)+go(a 2 +c u b 2 +db(l+c u )) 
1+90 (a+b) 



B. DeMarco, S. B. Papp, and D. S. Jin, Phys. Rev. 
Lett. 86, 5409 (2001); B. DeMarco and D. S. Jin, Sci- 
ence 285 1703 (1999); A. G. Truscott, K. E. Strecker, 
W. I. McAlexander, G. B. Partridge and R. G. Hulet, 
Science 291, 2570 (2001) ; F. Schreck. L. Khaykovich, 
K. L. Corwin, G. Ferrari, T. Bourdel, J. Cubizolles 
and C. Salomon, Phys. Rev. Lett. 87, 080403 (2001); 
Z. Hadzibabic, C. A. Stan, K. Dieckmann, S. Gupta, 



M. W. Zwier lein, A. Gorlitz and W. Ketterle, cond 



mat/01 12425 



S. R. Granade, M. E. Gehm, K. M. O'Hara 



[2] 



and J. E. Thomas, cond-mat/0111344 



[3] 

[4] 
[5] 

[6] 

[7] 

[8] 

[9] 

[10] 
[11] 

[12] 

[13] 

[14] 
[15] 



H. T. C. Stoof, M. Houbiers, C. A. Sackett, and R. G. 
Hulet, Phys. Rev. Lett. 76, 10 (1996); M. Holland, S. J. J. 
M. F. Kokkelmans, M. L. Chiofalo, and R. Walser, Phys. 
Rev. Lett. 87, 120406 (2001); L. Viverit, S. Giorgini, L.P. 
Pitaevskii, and S. Stringari, Phys. Rev. A 63, 033603 
(2001); H. Heiselberg, C. J. Pethick, H. Smith and L. 
Viverit, Phys. Rev. Lett. 85, 2418 (2000). 
P. Torma, and P. Zoller, Phys. Rev. Lett. 85, 487 (2000); 
G. M. Bruun, P. Torma, M. Rodriguez, and P. Zoller, 
Phys. Rev. A 64, 033609 (2001). 

Gh.-S. Paraoanu, M. Rodriguez and P. Torma, J. Phys. 
B 34, 4763 (2001). 

W. Zhang C. A. Sackett, and R. G. Hulet, Phys. Rev. A 
60, 504 (1999); J. Ruostekoski, Phys. Rev. A 60, R1775 
(1999), J. Ruostekoski, Phys. Rev. A 61, 033605 (2000); 

F. Weig and W. Zwerger, Europhys. Lett. 49, 282 (2000). 
M. A. Baranov and D. S. Petrov, Phys. Rev. A 62, 
041601 (R) (2000); M. Farine et al, Phys. Rev. A 62, 
013608 (2000); G. M. Bruun and C. W. Clark, J. Phys. 
B 33, 3953 (2000). 

D. M. Stamper-Kurn, A. P. Chikkatur, A. Gorlitz, S. In- 
noye, S. Gupta, D. E. Pritchard, and W. Ketterle, Phys. 
Rev. Lett. 83, 2876 (1999). 

D. Pines and P. Nozieres, The Theory of Quantum Liq- 
uids: Vol I (Benjamin, New York, 1966) 

A. Minguzzi, G. Ferrari and Y. Castin, Eur. Phys. J. D, 
17, 49 (2001). 

R. Cote and A. Griffin, Phys. Rev. B 48, 10404 (1993). 
P. W. Anderson, Phys. Rev. 112, 1900 (1958); N. N. Bo- 
goliubov, V. V. Tolmachev and D. V. Shirkov, A New 
Method in the Theory of Superconductivity (Academy of 
Science, Moscow 1958, New York 1959). 

B. W. Shore, The theory of coherent atomic excitation 
(Wiley, New York, 1990). Vols. I and II. 

L. P. Kadanoff and G. Baym, Quantum Statistical Me- 
chanics (Benjamin, New York, 1962). 
A. J. Legget, Phys. Rev. 147, 119 (1966). 

G. M. Bruun and B. R. Mottelson, Phys. Rev. Lett. 87, 
270403 (2001). 



G 



